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ABSTRACT 

We introduce a data reduction package written in Interactive Data Language (IDL) for the Magellan 
Echellete Spectrograph (MAGE). MAGE is a medium-resolution {R ~ 4100), cross-dispersed, optical 
spectrograph, with coverage from ~ 3000 - 10000 A. The MAGE Spectral Extractor (MASE) incor- 
porates the entire image reduction and calibration process, including bias subtraction, flat fielding, 
wavelength calibration, sky subtraction, object extraction and fiux calibration of point sources. We 
include examples of the user interface and reduced spectra. We show that the wavelength calibration 
is sufficient to achieve ~ 5 km s~^ RMS accuracy and relative fiux calibrations better than 10%. A 
light-weight version of the full reduction pipeline has been included for real-time source extraction 
and signal-to-noise estimation at the telescope. 

Subject headings: stars: low mass — stars: fundamental parameters — stars: M dwarfs — stars: 
luminosity function — stars: mass function — Galaxy: structure 
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1. INTRODUCTION 

Cross-dispersed spectrographs are attractive instru- 
ments for a variety of science applications, primarily for 
their broad wavelength coverage and high spectral res- 
olution. While these instruments maximize wavelength 
coverage and spectral resolution, the reduction and cali- 
bration of their data is complex. The spatial and disper- 
sion axes are typically not aligned with the CCD columns 
and rows; a given wavelength can span several columns 
on the chip. Furthermore, the curvature and tilting of 
the orders complicates the tracing of order edges and the 
extraction of spectra. 

Fortunately, many of these issues have been addressed 
by modern reduction techniques (e.g., Moreno et al. 1982; 
Rossi et al. 1985; Ponz et al. 1986; Goodrich & Veilleux 
1988; Marsh 1989; Mukai 1990; Hah et al. 1994; Piskunov 
& Valenti 2002; Kelson 2003). Pipelines exist and are 
regularly employed for extracting cross-dispersed opti- 
cal (Piskunov & Valenti 2002) and infrared (Gushing 
et al. 2004) data. Large, multi-fiber spectroscopic sur- 
veys, such as the Sloan Digital Sky Survey (SDSS; York 
et al. 2000; Stoughton et al. 2002; SubbaRao et al. 2002), 
the Large Sky Area Multi-Object Fiber Spectroscopic 
Telescope (LAMOST; Cui et al. 2008) and the Radial 
Velocity Experiment (RAVE; Steinmetz et al. 2006) also 
use pipeline reduction techniques that are directly appli- 
cable to the analysis of cross-dispersed data. 

We have developed a new reduction package for the 
Magellan Echellete spectrograph (MAGE; Marshall et al. 
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2008). The MAGE Spectral Extractor (MASE) is writ- 
ten in the Interactive Data Language (IDL) and is based 
on the SDSS spectro2d pipeline^, the Magellan Inamori 
Kyocera Echelle (MIKE; Bernstein et al. 2003) IDL 
pipeline^ (Bernstein, Buries & Prochaska, in prep) and 
the XIDL astronomy package^. The software employs 
a graphical user interface (GUI) built with IDL widgets 
to control the reduction process. A faster, "quicklook" 
version of the pipeline is packaged with MASE, for use 
at the telescope. This version can be used for on-the- 
fiy reductions and provides estimates of signal-to-noise 
ratios and spectral morphology. 

In this paper, we describe the details of the MASE 
reduction pipeline and GUI. In §2, we review the proper- 
ties of the MAGE spectrograph. The reduction pipeline 
is detailed in §3, with examples of the MASE GUI and 
reduced spectra. Our conclusions are reported in §4. 

2. MAGE 

The MAGE spectrograph is described in detail by Mar- 
shall et al. (2008), but we briefiy summarize details criti- 
cal to the reduction and calibration process here. MAGE 
is a medium-resolution {R :^ 4100, 22 km s~^ pixel"^) 
optical spectrograph currently mounted on the central 
folded port of the Clay telescope at the Las Campanas 
Observatory in Chile. Similar to the Echellette Spectro- 
graph and Imager (ESI; Sheinis et al. 2002) on the Keck 
II telescope, MAGE employs a grating in combination 
with two prisms to provide cross dispersion. The prisms 
are aligned to minimize anamorphic magnification. The 
instrument's wavelength coverage spans 3000 - 10500 A 
over 15 orders (6 < m < 20), continuous except for a 15 
A gap near 9500 A between orders six and seven. Five 
non-overlapping slits are available, with widths from 5.0'' 
to 0.5'', corresponding to resolutions of R ^ 1000 (~ 120 
km s-^) to R ^ 8000 (- 12 km s-^). Each slit is 10" 
long. The slit is imaged with a //1. 4 Schmidt vacuum 

^ Available at http://spectro.princeton.edu/ 

^ Available at http://web.mit.edu/~burles/www/MIKE/ 

^ Available at http://www.ucolick.org/~xavier/IDL/index.html 
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camera equipped with an E2V back-illuminated CCD. 
The chip is composed of 2048 x 1024 13.5 /im pixels. 
Each pixel subtends 0.33'' on the sky. Dispersion on the 
chip ranges from 0.3 A pixel" ^ (blue) to 0.6 A pixel" 
(red). A typical MAGE science observation is shown in 
Figure 1. 

Two calibration sources are incorporated into MAGE. 
A Thorium- Argon cathode tube is used for wavelength 
calibration, producing arc lines over the entire spectral 
window. ThAr arcs should be acquired after each science 
target to establish an accurate wavelength solution. A 
Xenon-flash bulb is used for fiat fielding in the bluer or- 
ders (m > 14), but this lamp cannot be used to fiat field 
the redder orders due to the presence of strong emission 
lines. We recommend the following calibration images 
to construct pixel and illumination fiats. Pixel fiats are 
constructed from three sets of observations. First, for 
the bluest orders (18 - 20), in focus twilight fiats should 
be acquired with the 5.0'' slit. This wide slit diffuses 
the Fraunhoffer lines in the sky spectrum. For the green 
orders (14 - 17), the Xe-fiash lamp observed with the sci- 
ence slit provides suitable illumination for fiat fielding. 
Finally, in the red orders (6 - 13), the bright quartz lamp 
and dome screen observed with the science slit provide 
suitable illumination of the optical path. For each set of 
orders, at least 10 focused fiat-field images are recom- 
mended to obtain a reliable median. Illumination fiats 
should also be acquired during twilight with the science 
slit. Since the prisms and gratings within MAGE do not 
move, the instrument is very stable and only one fiat 
field is typically needed per setup per night. If slits are 
changed during the night, a new set of fiats should be 
acquired after the change. Flux calibration is obtained 
by observing spectrophotometric standards (e.g., Oke 
1990; Hamuy et al. 1992, 1994; Massey et al. 1988). At 
least one standard per setup should be observed during 
the run. Note that observations of fiux calibrators are 
also required for object tracing. A table describing the 
necessary MAGE calibrators is listed below (Table 1). 



3. DATA REDUCTION 

3.1. User Interface 

The MASE GUI is based on IDL widgets. It is com- 
posed of two parts, the message window, which provides a 
log of the user's file selections and procedures executed, 
and the tabs, which allow for easy navigation through 
each step of the reduction process. An example of the 
MASE GUI is shown in Figure 2. Within the tabs, the 
user specifies the current working directory, constructs 
calibration images, and executes the main reduction soft- 
ware. Most common options, such as choosing which 
target to reduce, are available within the GUI interface. 
Other options are accessible only from the underlying 
IDL code, which can be executed from the command line. 
The default choices are optimized for point source extrac- 
tion. The MASE GUI and associated IDL programs are 
publicly available^^. 

10 The MASE website is located at 
htt p : / / web . mit . edu /j j b/ www / mase . html 



3.2. Organization and Structure 

For each night of observations, MASE generates an 
IDL data structure to organize the reduction. When the 
structure is first created, the program attempts to guess 
each type of exposure from header information and de- 
rived quantities. Each science target is assigned an ob- 
ject ID number and associated with the nearest temporal 
ThAr arc. The structure is written out to a FITS file, 
and can be inspected and edited with utilities included 
within MASE. 

The user selects the directory containing the raw 
MAGE frames to begin the reduction process. This selec- 
tion is echoed in the MASE message window. This data 
directory should be placed within a working directory 
(i.e. for each night's observing). A series of subfolders 
are created in directories parallel to the raw directory to 
house files generated during the reduction process. 

3.3. Calibration: Bias-subtraction 

The data reduction pipeline, outlined in Figure 3, be- 
gins with bias subtraction of every raw frame. The me- 
dian bias level is measured from a 128 x 128 pixel over- 
scan region on the MAGE detector and subtracted from 
the entire image. Each image is then trimmed to 2048 x 
1024 pixels and pixel values are converted from ADUs to 
electrons, using a previous empirical measurement of the 
gain, which is recorded in the image header. The prede- 
termined read noise and electron counts (shot noise) are 
used to construct an inverse variance image, which is a 
first estimate of the uncertainty in the counts per pixel 
position caused by photon counting noise and detector 
read noise. 

3.4. Calibration: Flat fielding and Order Definition 

A fiat-field image is constructed for each illumination 
source described in §2. First, a median box filter is ap- 
plied to each fiat, smoothing the image over 30 pixel 
scales and removing large scale variations. Each un- 
smoothed fiat field image is then normalized by its me- 
dian filtered counterpart. The normalized fiat-field im- 
ages are then coadded and further filtered, where outly- 
ing pixels with deviations larger than two standard devi- 
ations from the mean are not included in the final coadd. 
Next, the separate red, green and blue fiats are merged 
into a single "super fiat" . The red fiat is used for orders 
6-13, the green fiat spans orders 14-17, and the blue fiat 
covers orders 18 through 20. For typical cross-dispersed 
fiat fields, the order edges can exhibit sharp discontinu- 
ities. Each order was trimmed by four pixels on its upper 
and lower bound to minimize this effect. Furthermore, 
inter-order pixels were set to unity. An example of a su- 
per fiat is shown in Figure 4. The imaged is stretched to 
greatly enhance the contrast of the image and highlight 
features within the fiat-field image. In orders 6 and 7, 
fringing at the level of ^ 10% is clearly visible. Orders 
8 through 14 are very fiat, showing no almost no vari- 
ations, except for a few small "holes" (i.e., such as the 
one visible in order 12). Orders 16 and 17 display some 
variations due to small fluctuations in the anti-refiective 
coating on the CCD. Finally, order 20 is set to unity, as 
obtaining the required number of counts (~ 2-3 x 10^) at 
these wavelengths (< 3300 A) for fiat fielding is severely 
hampered by atmospheric extinction. The illumination 
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Fig. 1. — An example MAGE spectrum of a flux standard. The central wavelengths (blue) and order numbers (red) are listed near each 
order. The position of the wavelength labels do not correspond to their actual location on the chip. Note that the sensitivity of MAGE 
falls off sharply in orders 6 and 20. 
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Raw path has been updated to /Volu[iies/Beefnut/v2Reductions/0ct2008/ut081006/Raw/ 
Current working directory is/Volurnes/Beefnut/v2Reductions/0ct2008/ut0S1006 
Mage structure created and written to mage *f its 
Mage structure has been edited and written to mage+fits 

Green flat files added as /Volumes/Beef nut/v2Reductions/0ct2008/ut0S1006/Raw/mage0021, fits 
Green flat files added as /Volumes/Beef nut/v2Reductions/0ct2008/ut08100G/Raw/mage0022*f its 
Green flat files added as /Volumes/Beef nut/v2Reductions/0ct2008/ut081006/Raw/mage0023, fits 
Green flat files added as /Volumes/Beef nut/v2Reductions/0ct2008/ut081006/Raw/mage0024. fits 
Green flat files added as /Volumes/Beef nut/v2Reductions/0ct2008/ut08100G/Raw/mage0025, fits 
Green flat files added as /Volumes/Beefnut/v2Reductions/0ct2008/ut081006/Raw/mage0026,f its 
Green flat files added as /Volumes/Beef nut/v2Reductions/0ct2008/ut08100G/Raw/mage0027*f its 
Green flat files added as /Volumes/Beef nut/v2Reductions/0ct2008/ut081006/Raw/mage0028, fits 
Green flat files added as /Volumes/Beef nut/v2Reductions/0ct2008/ut081006/Raw/mage0029. fits 
Green flat files added as /Volumes/Beef nut/v2Reductions/0ct2008/ut08100G/Raw/mage0030, fits 
Red flat files added as /Volumes/Beef nut/v2Reductions/0ct2008/ut081006/Raw/mage0071, fits 
Red flat files added as /Volumes/Beef nut/v2Reductions/0ct2008/ut08100G/Raw/mage0072. fits 
Red flat files added as /Volumes/Beef nut/v2Reductions/0ct2008/ut081006/Raw/mage0073, fits 
Red flat files added as /Volumes/Beef nut/v2Reductions/0ct2008/ut081006/Raw/mage0074. fits 
Red flat files added as /Volumes/Beef nut/v2Reductions/0ct2008/ut08100G/Raw/mage0075. fits. 
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Fig. 2. — The MASE GUI. It is composed of the message window and tabs. The message window, located in the upper part of the GUI, 
provides feedback and conflrmation of user actions. Underneath the message window, the tabs separate the major steps in the calibration 
process and allow the user to select which flies are calibrated and reduced. The Flats tab is shown in this example. 
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TABLE 1 
Suggested Calibrations 



Calibration Source 


Slit 




Applicable Orders 


Exposure Time^ 


Num. of Exp. 


Quartz Lamp 


Science 


Slit 


6-13 


40s 


>io 


Xe-Flash Lamp 


Science 


Slit 


14-17 


35s 




Twilight 






18-20 


- 60s b 




Twilight^ 


Science 


Slit 


6-20 


- 60s 




Flux Calibrator 


Science 


Slit 


6-20 


Varies 




Thorium Argon 


Science 


Slit 


6-20 


3s 


Varies 



^Exposure times are given for the 0.7^^ and should be scaled down for larger slits. 
^Twilight exposure times should be adjusted for changing sky conditions. 
*^ These exposures are used for illumination flats. 



Raw MAGE data 



Arcs 



Extract Arcs 



Locate and fit 
lines 



Fit tilt in lines 



X,Y Pixel to 
Wavelength 
map 



Overscan Subtract 



I- Flats 



Create Trace 



Create separate 
R,G,B flats 



Trace Orders 



Create 
superflat 



Vacuum & Heliocentric 
Corrections 



Flux Calibrate 



Combine Orders 



Final 1D Spectrum 



Science 



^ Divide by flat 



Subtract Sky 
(B-spline) 



Extract 
Spectrum 



Apply 
Wavelength map 



Fig. 3. — An illustrative flowchart of the MASE reduction pipeline. Each major step in the pipeline discussed in the text is shown. 
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TABLE 2 

Wavelength ranges by order 



Order 




^central 


^max 


g 


9500 


9893 


10285 


7 


8128 


8752 


9475 




7112 


7700 


8288 


9 


6320 


6844 


7368 


10 


5713 


6172 


6630 


11 


5225 


5626 


6027 


12 


4286 


4906 


5525 


13 


4469 


4785 


5100 


14 


4170 


4445 


4720 


15 


3905 


4148 


4390 


16 


3675 


3888 


4100 


17 


3467 


3660 


3853 


18 


3287 


3453 


3620 


19 


3187 


3308 


3430 


20 


3197 


3254 


3311 



correction is computed by fitting the variation in the spa- 
tial direction of a smoothed twihght flat fleld along the 
X coordinate of each order. 

A critical step in reducing cross-dispersed data is the 
identification and measurement of order edges. Trace 
flats are constructed by summing red, green and blue 
flat-fleld images for each slit used. The trace image is 
convolved with a sawtooth fllter, enhancing the order 
edges. A fourth order Legendre polynomial is then flt 
to the upper and lower edges of each order as a function 
of column (x pixel position). The flt coeflicients deflne 
the order edges, centers and widths on the chip. The 
order information is output to a FITS flle, and the flt is 
displayed for visual conflrmation by the user. 



3.5. Calibration: Wavelength Solutions 

A wavelength solution is computed from each ThAr 
arc image. An example arc frame is shown in Figure 
5. Because the dispersion of light is not aligned along 
rows and columns, each pixel corresponds to a speciflc 
wavelength. The wavelength solution is derived in two 
parts. First, the pixels at the center of each order are 
mapped to wavelength. A boxcar extraction down the 
center of each order is used to extract a ID arc spectrum 
in air wavelengths, which is compared to a ThAr line at- 
las (Murphy et al. 2007). We note that many currently 
available line lists are appropriate for higher resolution 
echelle spectrographs (i.e., MIKE). To construct a clean 
line list for MAGE, we removed line blends by eye, us- 
ing the line atlas of Murphy et al. (2007) smoothed to 
MAGE resolutions as a guide. An archived solution was 
produced by manually identifying emission lines in every 
order. A Fourier transform is used to calculate a shift be- 
tween the extracted arc spectrum and the archived wave- 
length solution. These shifts are incorporated within an 
initial guess at the wavelength solution. The line cen- 
troids of the extracted arc spectrum are identifled and 
a low-order Chebyshev polynomial is flt to A vs. col- 
umn {x pixels) using least squares. This result is used 
to identify slightly weaker lines excluded in the original 
flt, and the polynomial is re-flt. Next, a two-dimensional 
solution is computed, fltting the product of the wave- 
length and echelle order as a function of CCD row and 



order number. The 2D solution allows for interpolation 
through orders with a small amount of clean ThAr lines, 
such as orders 6 and 20. It also constrains the wavelength 
solution to be the same in regions of overlap, where the 
ID solutions are totally independent. The 2D solution is 
not complete, as it only maps the center of each order to 
wavelength. The tilt of the emission lines, which varies 
with position on the chip, must be incorporated into the 
flnal wavelength calibration. In each order, strong emis- 
sion lines are flt with a straight line, and the resulting 
slope is used to determine the tilt of the arcs as a func- 
tion of position. This information is combined with the 
center-order wavelength solution, and each science pixel 
is mapped to wavelength. A 2D FITS image of the pixel- 
wavelength map is written in the Arcs directory. 

Multiple quality assurance (QA) plots are produced 
during the wavelength calibration process. This allows 
the user to inspect both the one dimensional wavelength 
solution for each order and the two dimensional flts. The 
root mean square (RMS) deviation in pixels is recorded 
for each order, and lines that are rejected from the flnal 
flt are flagged. Example QA plots are shown in Figures 
6 and 7. Typical RMS values for the one dimensional 
flts are ~ 0.05, and ranges up to 0.2 pixels for orders 
with very few ThAr lines (6 & 20). This corresponds to 
roughly 0.02 - 0.07 A and sets the velocity precision of the 
MAGE spectrograph. This precision is approximately 
5 km s~^, as verifled by cross-correlating observations 
of known TW Hydrae Association members using the 
0.7'' slit with known radial velocity standards (Bochanski 
et al. 2007). 

MAGE is an exceptionally stable instrument, espe- 
cially over the course of a night of observations. To test 
of the wavelength solution stability, we extracted 13 sky 
spectra from a given night of observations. No heliocen- 
tric corrections were applied, and the telluric emission 
lines were cross-corellated against each other. The rel- 
ative shifts over course of the night were typically less 
than 1 km s~^, with a mean of > 0.3 km s~^. This shift 
corresponds to > 0.01 pixel, indicating MAGE is a very 
stable instrument. Thus, the observer can confldently 
apply an arc taken before or after an observation, but 
can avoid the overhead cost of two arcs. 

3.6. Extraction: Sky-Subtraction, Object Finding, 
Aperture Tracing and Extraction 

Using the calibration images described above, each 
science target is overscan subtracted and flat flelded. 
Yet, night sky (telluric) emission lines and background 
remain. Most traditional sky-subtraction algorithms 
model telluric features after rectifying the observed spec- 
tra onto a uniform wavelength grid. Unfortunately, re- 
binning the spectrum correlates the noise in the im- 
age, combining pixels that were previously independent. 
However, recent data reduction pipelines (e.g., Bolton & 
Buries 2007) have introduced basis splines (B-Splines; 
Kelson 2003), or piecewise polynomial solutions deflned 
by an order and breakpoints to produce a model of the 
sky emission and remove it from the science image. At 
each breakpoint, the B-spline is continuous through n — 2 
derivatives, where n is the order. Typically, 4th order B- 
splines are employed for sky subtraction, which use cubic 
piecewise functions. B-splines are flexible enough to flt 
the night sky emission features as a function of CCD 
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3453 



3254 



20 



Fig. 4. — An example "superflat" combining separate red, green and blue flats. The image has been stretched to enhance the contrast. 
As in Figure 1 the order numbers are given in red and the central wavelengths of each order are listed in blue. There is signiflcant fringing 
in the reddest orders, and note that order 20 has all pixels set to unity. 
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Fig. 5. — Shown is an example ThAr arc spectrum. Note the dearth of emission lines in the reddest and bluest orders. 



pixels, rather than wavelength. This greatly simplifies 
the error analysis and propagation, and results in near 
Poisson-limit sky subtraction. Additional discrepant fea- 
tures, such as cosmic rays or bad pixels are easily iden- 
tified and removed from the model sky image. 

To start the extraction process, the object is crudely 
traced and masked. Each order is collapsed along wave- 
length, and the spatial center of the profile is located with 
a median. The region along this crude trace is masked 
out, and a sky model is constructed from the background. 
For each science target, a 4th order B-spline is fit as func- 
tion of pixel position along a trace offset from the order 
center. Using the tilt mapped during the wavelength cal- 
ibration, the initial sky model is propagated across the 
entire spatial extent of the slit. This initial sky model 
is used as a best guess for the final sky, which is later 



output as a FITS image in the Final directory. 

Next, an object finding code is run to refine the crude 
object trace above. As described above, the code uses 
the trace of a previously reduced flux standard as a best 
guess for the object's position. Since the flux standard is 
bright, locating and extracting its spectrum is straight- 
forward and does not require a precise initial guess of its 
position on the chip. First, the trace of each object is 
compared to the standard star trace in each order, and 
the distance between the science target and flux standard 
in each order is computed. The traces are grouped by a 
friends-of-friends algorithm. For orders in which a trace 
was not detected (i.e., a Lyman-limit system or T dwarf) 
the median distance from the standard star trace (for all 
detected orders) is used. Note that this does not account 
for any shift due to atmospheric refraction between the 
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Fig. 6. — An example of the QA plots produced during ID wavelength calibration. Plotted in each panel are the residuals (in pixels) 
as a function of wavelength, for each order. The black crosses correspond to ThAr emission lines that were used in the final fit, and red 
asterisks are lines that were rejected. The order number, RMS (in pixels), dispersion (A pixel" ■'^) and polynomial order is reported for each 
panel. Note that the typical RMS is less than 0.1 pixels and usually ~ 0.05 pixels. 
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Fig. 7. — An example of the QA plots produced during 2D wavelength calibration. The blue line is the 2D Chebyshev polynomial fit, 
and the crosses are the individual line centers, shown as a function of column (x pixel position) and wavelength. This step improves the 
accuracy in orders with few ThAr lines, such as order 6 (the rightmost order). 



standard and science object. Once the trace center has 
been determined, the flux-weighted mean position as a 
function of CCD row is measured, and a low-order Leg- 
endre polynomial is fit to trace on the chip. A structure, 
similar to the one describing the order edges, is generated 
to describe the location of the object along the chip. 

After finding the object, the dispersion axis is collapsed 
and the counts are plotted against fractional slit height. 
An example of the collapsed slit profile is shown in Fig- 
ure 8. As the imaged slit width in pixels varies between 
orders, identifying the object in relative slit coordinates 
simplifies the extraction code. Once the object trace is 
located on the chip, a two-step process is used to mea- 
sure its flux. First, a simple boxcar extraction is used 
to quickly estimate the signal-to-noise ratio (SNR) in 



each order. The orders are then ranked by SNR, and 
an optimal extraction is used to extract the flux from 
the point source spectrum. Optimal extraction was in- 
troduced by Horne (1986) and Robertson (1986) and up- 
dated for use in cross dispersed data by Mukai (1990) and 
Marsh (1989). Briefly, optimal extraction weights each 
pixel by the fraction of the total flux impinging on that 
pixel. These weights are calculated from the slit profile, 
which is measured along the collapsed order (seen in Fig- 
ure 8). The profile is fit with a B-spline for high SNR 
observations, and for low SNR orders (with SNR < 2.49), 
the profile parameters determined in the higher SNR or- 
ders are assumed. If the median SNR of the spectrum 
falls below the limit, a Gaussian is assumed. Note that 
optimal extraction is ideal for spectroscopy of unresolved 
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point sources, and will recover the proper profile for ex- 
tended objects as well, provided their slit profile shape 
does not vary strongly with wavelength. 

Extraction takes place after the object profile and op- 
timal extraction weights have been defined. B-splines 
are employed to simultaneously model both the sky and 
object spectra. This routine improves the model sky es- 
timate described above. After the first B-spline fit, pixels 
with a deviation greater than 3.5 sigma from the model 
are identified. This mask is designed to remove bad pix- 
els and cosmic ray hits. The fit is repeated, with out- 
liers masked out. The two dimensional sky and object 
images are recorded, and the spectrum of the object is 
written out to a file. The 2D images can be inpected 
with a separate IDL routine (mageJook) , which overlays 
image masks, order edges and the object trace on the 
2D image. The wavelength solution determined above 
is adopted during the extraction, after applying a helio- 
centric correction and converting to vacuum coordinates. 

3.7. Post Extraction: Combining Orders and Flux 
Calibration 

After extraction, each order's spectrum is recorded in 
counts vs. wavelength. The spectra are then rebinned 
onto a common logarithmic wavelength grid. Prior to 
this step, each pixel on the chip has been treated inde- 
pendently, which avoids correlated errors. However, a 
common wavelength grid is necessary for co- adding mul- 
tiple exposures, which occurs after rebinning. Diff"erent 
exposures of a common science target are co-added by 
default, but the user has the option of recording individ- 
ual spectra, useful for monitoring a given target (i.e., for 
radial velocity variations). 

After co-addition, a fiux calibration is applied to the 
data. The user can specify an archival solution or de- 
termine a new sensitivity function. The standard star 
spectrum is extracted as described above, then compared 
to spectrophotometric atlases incorporated with MASE 
(Oke 1990; Hamuy et al. 1992, 1994; Massey et al. 1988). 
A list of the standards included can be found on the 
MASE website. The user can interactively mask out bad 
regions or strong absorption lines in the flux calibrator's 
spectrum. A sensitivity function is determined (counts 
/ fiux as a function of wavelength) and applied to the 
science spectrum. Typical relative fluxes are good to ^ 
10%, as shown in Figure 9. A FITS flle containing the 
spectrum for each order is output to the Object directory. 

Finally, each co-added spectrum is combined to a one- 
dimensional fiux vs. wavelength spectrum. At the over- 
lapping edges of sequential orders, the average is com- 
puted and recorded for that wavelength bin. Agreement 
between the orders is usually good to < 5%, as shown in 



Figure 10. The ID flux and error spectra are output to 
the FSpec directory. 

3.8. Quicklook Tool 

A "slimmed-down" version of the full pipeline is in- 
cluded with MASE for use at the telescope. This quick- 
look tool (QT) can process a raw CCD image in a 
few minutes and allow the observer to quickly estimate 
signal-to-noise and inspect spectral morphology. It em- 
ploys archival order edge traces, wavelength solutions 
and sensitivity functions. The QT uses boxcar extrac- 
tion instead of optimal extraction, resulting in major im- 
provements in speed. However, emission lines are com- 
monly mis-identified as night sky lines and subtracted 
from the spectrum. The user should always use the full 
reduction package for scientific analysis. The QT is ac- 
cessed by clicking on the "Quicklook" tab in the MASE 
GUI. The only user input is the filename of the science 
exposure. If the user does not have MASE installed on 
their own computer, the QT is also installed locally at 
Las Campanas. 

4. SUMMARY 

We have developed a new data reduction package 
for the MAGE spectrograph. The MASE GUI and 
pipeline is written in IDL and incorporates modern sky- 
subtraction and optimal extraction techniques. The pro- 
gram is available for download online. The pipeline is de- 
signed for reductions of unresolved point sources. How- 
ever, for spectra of extended objects, MASE can still be 
used to bias-subtract, fiat field and wavelength calibrate 
the observations. Future versions of MASE will include a 
telluric correction routine, which will require additional 
calibration observations, and scripting capabilities for re- 
ducing multiple nights of data. 
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